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We study the nearest neighbor XXZ Heisenberg quantum antiferromagnet on the kagome lattice. 

Here we consider the effects of several perturbations: a) a chirality term, b) a Dzyaloshinski-Moriya 
term, and c) a ring-exchange type term on the bowties of the kagome lattice, and inquire if they 
can support chiral spin liquids as ground states. The method used to study these Hamiltonians 
is a flux attachment transformation that maps the spins on the lattice to fermions coupled to a 
Chern-Simons gauge field on the kagome lattice. This transformation requires us to consistently 
define a Chern-Simons term on the kagome lattice. We find that the chirality term leads to a chiral 
spin liquid even in the absence of an uniform magnetic field, with an effective spin Hall conductance 
of a’fy = I in the regime of XY anisotropy. The Dzyaloshinski-Moriya term also leads a similar 
chiral spin liquid but only when this term is not too strong. An external magnetic field also has 


the possibility of giving rise to additional plateaus 
XY regime. Finally, we consider the effects of a 
coupling constant is large enough, it may trigger a 
spontaneous breaking of time-reversal invariance. 

I. INTRODUCTION 

Quantum Heisenberg antiferromagnets on two- 
dimensional kagome lattices are strongly frustrated 
quantum spin systems making them an ideal candidate 
to look for exotic spin liquid type states. This model has 
been the focus of many theoretical and numerical efforts 
for quite some time and, in spite of these efforts, many 
of its main properties remain only poorly understood. 

On a parallel track, there has also been significant 
progress on the experimental side with the discovery of 
materials such as Volborthite and Herbertsmithite whose 
structures are closely represented by the Heisenberg 
antiferromagnet on the kagome lattice. An extensive 
analysis of the present experimental and theoretical 
status of this problem can be found in a recent review 
by Balents.^ 

It is well known that frustrated antiferromagnets can 
give rise to magnetization plateaus in the presence of an 
external magnetic field. Although the nature of these 
plateaus depend on the type of model as well as the 
regime under study, in two dimensions these states are 
expected to be topological phases of the chiral spin liquid 
type in the XY regime.!^ Previous works have also shown 
that in the Ising regime these magnetizat ion plateaus 
behave like valence bond crystals (VBC).I^ In a re¬ 
cent work, the same authors studied the nearest neigh¬ 
bor XXZ Heisenberg model on the kagome lattice using 
a newly developed flux attachment transformation that 
maps the spins which are hard-core bosons to fermions 
coupled to Chern-Simons gauge field.l^ Here we showed 
that in the presence of an external magnetic field, this 
model gives rise to magnetization plateaus at magnetiza- 


which also behave like chiral spin liquids in the 
ring-exchange term and find that, provided its 
phase transition into a chiral spin liquid by the 


tions m = ^, I and | in the XY regime. The plateaus 
at m = I and m = | behave like a Laughlin fractional 
quantum Hall state of bosons with an effective spin Hall 
conductance of a^y = | whereas the plateau at | was 
equivalent to a Jain state with a^y = |. 

In spite of the considerable work done on this prob¬ 
lem since the early 1990s, the kagome Heisenberg an¬ 
tiferromagnet without an external magnetic field is less 
well understood. Density matrix renormalization group 
(DMRG) studies (which crucially also included sec¬ 
ond and third neighbor neighbor antiferromangnetic ex¬ 
change interactions) have found a topological phase in 
the universality class of the Z 2 spin liquid.l^HSl A Z 2 spin 
liquid has been proposed by Wang and VishwanatlP (us¬ 
ing slave boson methods), and by Fisher, Balents and 
Girviiiinini, 

in generalized ferromagnetic XY model with 
ring-exchange interactions. Similarly, a Z 2 spin liquid 
was shown to be the ground state for the kagome anti¬ 
ferromagnet in the quantum dimer approximation.^^ We 
should note, however, that an entanglement renormaliza¬ 
tion group calculatioiP^ appears to favor a complex 36 
site VBC as the ground state of the isotropic antiferro¬ 
magnet on the kagome lattice. 

A more complex phase diagram (which also includes a 
chiral spin liquid phase as well as VBC, as well as Neel 
and non-colinear antiferromagnetic phases) was subse¬ 
quently found in the same extension of the kagome anti¬ 
ferromagnet by Gong and coworkers.^ It has also been 
suggestecP^ that a model of the kagome antiferromagnet 
with second and third neighbor Ising interactions may 
harbor a chiral spin liquid phase. Still, variational wave 
function studies have a lso suggested a possible U{1) Dirac 
spin liquid state,but this is not consistent with the 
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DMRG results. Bieri and coworkerJi^ used variational 
wave function methods to study a kagome system with 
ferro and antiferromagnetic interactions, and find evi¬ 
dence for a gapless chiral spin liquid s tate. Other studies 
that used non-linear spin-wave theorjffsHni tiad predicted 
a quantum disordered phase for spin-1/2 kagome antifer- 
romagnets with XY anisotropy. 

The most recent DMRG studies on the plain (with¬ 
out further neighbors) kagome Heisenberg antiferromag- 
net give a strong indication of a gapped time-reversal 
invariant ground state in the Z 2 topological class.® As it 
turns out there are two possible candidates Z 2 spin liq¬ 
uids: the double Ghern-Simons theorjED with a diagonal 
AT-matrix, K = diag(2, —2), and the Kitaev toric cod^^ 
(which is equivalent to the deconfined phase of the Ising 
gauge theorjEH). An important recent result by Zaletel 
and VishwanatfP^l proves that the double Ghern-Simons 
gauge theory is not allowed for a system with translation 
and time-reversal invariance such as the kagome antifer- 
romagnet, which leaves the Z 2 gauge theory as the only 
viable candidate. 

In a recent recent and insightful work, Bauer and 
coworker^^ considered a kagome antiferromagnet with 
a term proportional to the chiral operator on each of 
the triangles of the kagome lattice. In this model time- 
reversal invariance is broken explicitly. These authors 
used a combination of DMRG numerical methods and 
analytic arguments to show that, at least if the time- 
reversal symmetry breaking term is strong enough, the 
ground state is a chiral spin liquid state in the same topo¬ 
logical class as the Laughlin state for bosons at filling 
fraction 1/2. This state was also found by u^in the 1/3 
and 2/3 magnetization plateaus. A more complex phase 
diagram was recently found in a model that also included 
second and third neighbor exchange interactions .ESI On a 
separate track, Nielsen, Cirac and SierraP^ constructed 
lattice models with long range interactions with ground 
states closely related to the 1 / = 1/2 Laughlin wave func¬ 
tion for bosons and, subsequently, deformed these models 
to systems with short range interactions on the square 
lattice that, in addition to first and second neighbor ex¬ 
change interactions, have chiral three-spin interactions 
on the elementary plaquettes, and showed (using finite- 
size diagonalizations on small systems) that the ground 
state is indeed a chiral spin liquid.ESl 

In this paper we investigate the occurrence of chiral 
spin liquid states in three extensions of the quantum 
Heisenberg antiferromagnet on the kagome lattice: a) by 
adding the chiral operator acting on the triangles of this 
lattice, b) by considering the effects of a Dzyaloshinski- 
Moriya interaction (both of which break time reversal 
invariance explicitly), and c) by adding a ring-exchange 
term on the bowties of the kagome lattice. This term, 
equivalent to the product of two chiral operators of the 
two triangles of the bowtie, does not break time-reversal 
invariance explicitly, and allows us to investigate the pos¬ 
sible spontaneous breaking of time-reversal. Through¬ 
out we will use a flux-attachment procedure suited for 


the kagome lattice that we developed recently.EEH This 
method is well suited to investigate chiral spin liquid 
phases but not as suitable for Z 2 phases (at least not 
straightforwardly). 

Throughout this paper we use a flux attachment trans¬ 
formation that maps the spins on the kagome lattice to 
fermions coupled to a Ghern-Simons gauge field.ESI The 
flux attachment transformation requires us to rigorously 
define a Ghern-Simons term on the lattice. Previously, 
such a lattice Ghern-Simons term ha d only been written 
down for the case of the square lattice.EiEII]\/[ore recently, 
we have shown how to write down a lattice version of the 
Ghern-Simons term for a large class of planar lattices.E^ 
Equipped with these new tools, we can now study the 
nearest neighbor Heisenberg Hamiltonian, chirality terms 
and the Dzyaloshinski-Moriya terms on the kagome lat¬ 
tice, as well as ring exchange terms. Furthermore, this 
procedure will allow us to go beyond the mean-field level 
and consider the effects of fluctuations which are gener¬ 
ally strong in frustrated systems. 

In spite of the successes of these methods, we should 
point out some of its present limitations. In a separate 
publication,E^ we showed that the discretized construc¬ 
tion of the Ghern-Simons gauge theory (on which we rely 
heavily) can only be done consistently on a class of planar 
lattices for which there is a one-to-one correspondence 
between the sites (vertices) of the lattice and its plaque¬ 
ttes. At present time this restriction does not allow us 
to use these methods to more general lattices of interest 
(e.g. triangular). 

We should also stress that in this approach there is no 
small parameter to control the accuracy of the approx¬ 
imations that are made. As it is well known from the 
history of the application of similar methods, e.g. to the 
fractional quantum Hall effects,E^lEZl tJ^at they can suc¬ 
cessfully predict the existence stable of chiral phases pro¬ 
vided the resulting state has a gap already at the mean 
field level. These methods predict correctly the univer¬ 
sal topological properties of these topological phases, in 
the form of effective low-energy actions, which encode 
the correct form of the topologically-protected responses 
as well as the large-scale entanglement properties of the 
topological phases.ESEH However, these methods cannot 
predict with significant accuracy the magnitude of di¬ 
mensionful quantities such as the size of gaps since there 
is no small parameter to control the expansion away from 
the mean field theory. Thus, to prove that the predicted 
phases do exist for a specific model is a motivation for 
further (numerical) studies. In particul^ the recent nu¬ 
merical results of Bauer and coworkerg^ are consistent 
with the results that we present in this paper. 

Using these methods, we find that chiral spin liquid 
phases in the same topological class as the Laughlin state 
for bosons at filling fraction 1/2 occur for both the model 
with the chiral operator on the triangles and for the 
Dzyaloshinski-Moriya interaction. We find that the chi¬ 
rality term opens up a gap in the spectrum and leads to a 
chiral spin liquid state with an effective spin Hall conduc- 
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tivity of a^y = ^ in the XY regime. This is equivalent 
to a Laughlin fractional quantum Hall state for bosons 
similar to the m = | magnetization plateau found in our 
earlier work.l^We also show that this chiral state survives 
for small values of the Dzyaloshinski-Moriya term. The 
main motivation for this study was motivated by a recent 
numerical worlP^ where the authors studied the nearest 
neighbor Heisenberg model in the presence of the chiral¬ 
ity and Dzyaloshinski-Moriya terms and found evidence 
of a similar fractional quantum Hall state with filling frac¬ 
tion i. The results we obtain agree qualitatively with 
the results obtained in the numerical work. Since this 
is the same state that we obtained at the 1/3 and 2/3 
magnetization plateaus, we searched for more complex 
topological phases by also adding a magnetic field. We 
also consider the effects of adding the chirality term and 
Dzyaloshinski-Moriya terms in the presence of an exter¬ 
nal magnetic field at the mean-field level. (For numerical 
studies of magnetization plateaus in the kagome anti- 
ferromagnets see the recent work by Capponi et al^. 
Once again, this is expected to give rise to magnetiza¬ 
tion plateaus. In the XY regime, we again find some 
of the same plateaus that were already obtained for just 
the case of the Heisenberg model.l^ In addition, we also 
find an additional plateau at m = ^ which has an ef¬ 
fective spin Hall conductivity of a^y = |. Finally, in 
order to address the possible spontaneous breaking of 
time-reversal, we consider the effects of a ring-exchange 
term on the bowties of the kagome lattice. Provided that 
the coupling constant of this ring-exchange term is large 
enough, this term triggers the spontaneous breaking of 
time-reversal symmetry and leads to a similar chiral spin 
liquid state with a^y = ^. 

The paper is organized as follows. In Sec. [Hj we 
will explicitly write down the XXZ Heisenberg Hamil¬ 
tonian and the chirality terms that we will consider. We 
will then briefly review the flux attachment transforma¬ 
tion that maps the hard-core bosons (spins) to fermions 
and write down the resultant Chern-Simons term on the 
kagome lattice. A more detailed and formal discussion 
on the flux attachment transformation and the issues 
related to defining a Chern-Simons term on the lattice 
is presented elsewhere.l^ An explicit derivation of the 
Chern-Simons term for the case of the kagome lat tice 
was also presented in our earlier work.l^ In Sec. IV 


we 


setup the mean-field expressions for the nearest neigh¬ 
bor Heisenberg model in the presence of a chirality term. 
We will then begin by analyzing the mean-field state in 
the absence of the chirality terms and analyze the states 
obtained in the XY and Ising regimes in Sec. |IV A| 
Then we will consider the effects of adding the chiral¬ 
ity term in Sec. |IVB[ At the mean-field level we will also 
study the effects of adding an external magnetic field to 
such a system by analyzing the Hofstadter spectrum in 
Sec. IV C[ Next, we will repeat the same analysis but 
with the Dzyaloshinski-Moriya term added to the XXZ 
Heisenberg model in Sec. |IVD[ In Sec. |Vj we will re¬ 
turn to the state discussed in Sec. IV B and expand the 


mean-field state around the Dirac points and write down 
a continuum version of the action. From here we will sys¬ 
tematically consider the effects of fluctuations and derive 
an effective continuum action. Finally, we will consider 
a model for spontaneous time-reversal symmetry break¬ 
ing by adding, to the XXZ Heisenberg antiferromagnetic 
Hamiltonian, a ring exchange term on the bowties of the 
kagome lattice in Sec |vg In Section m we summarize 
the key results obtained in the paper, and discusses sev¬ 
eral open questions. 


II. FLUX ATTACHMENT ON THE KAGOME 
LATTICE 


In this section, we will briefly review the theory of flux 
attachment on the kagome lattice in the context of the 
XXZ Heisenberg HamiltoniarP. This will set the stage 
for the use of the flux attachment transformation that we 
will use to study these models. 

We will begin with the nearest neighbor XXZ Heisen¬ 
berg Hamiltonian in the presence of an external magnetic 
field 


H = J ^ {5/5; + 5/5J + A5/5|} - ^ 5/ 

(*J> i 

( 2 . 1 ) 

where (i,j) refer to the nearest neighbor sites on the 
kagome lattice and A is the anisotropy parameter along 
the z direction, refers to the external magnetic 

field. Using the flux attachment transformation, the 
spins (which are hard-core bosons) can be mapped to 
a problem of fermions coupled to a Chern-Simons gauge 
field. The resultant action takes the form 


5 = 5f(^,^*, A^) + 5„t(A^) + 05 cs(A^) (2.2) 


The 5^5^ and 5^5^ terms map to the fermionic hop¬ 
ping part in the presence of the Chern-Simons gauge field, 
Aj(x), and 5*5^ terms map to fermionic interaction term 
as shown in the below equations 


SF{^,r,A^) = /?[ tp*{x) {iDo +Ai)V'(a 


{x,x') 



n{x') 


(2.3) 

where Dq = Sq -I- zAg is the covariant time derivative, 
(x, x') stands for nearest neighbor sites x and x' on the 
kagome lattice and the space-time coordinate x = (x, t). 
The temporal gauge fields Ag live on the sites of the 
kagome lattices and the spatial gauge fields A^ (x) live on 
the links of the lattice as can be seen in the unit cell of 
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FIG. 1: Unit cell of the kagome lattice. The unit cell 
has three sub-lattice sites (labelled a, b and c) and 
three plaquettes (two triangles and one hexagon). The 
flux attachment transformation proceeds by attaching 
the fluxes in each of the plaquettes to its corresponding 

sites. 


the kagome lattice in Fig[T] 

The density operator n{x) = ip*{x)tp{x) is related to 
the spin component as follows 

= ^ - n{x) (2.4) 

The above expression also allows us to absorb the ex¬ 
ternal magnetic held term [hext] in to the dehnition of 
the chemical potential /i, i.e. in the fermionic language 
the effect of the external magnetic held can be mimicked 
by changing the fermion density on the lattice. For a 
majority of this paper, we will focus on the case where 
hext = 0. This would correspond to the case of half- 
hlling in the fermionic theory after the hux attachment 
transformation. 

Now all that remains is the Chern-Simons term on the 
kagome lattice. An explicit derivation of this term for 
the case of the kagome lattice was already presented in 
an earlier paper.^ A more detailed and rigorous repre¬ 
sentation of a Chern-Simons term on generic planar lat¬ 


tices is also presented elsewhere.l^ Here, we will simply 
reproduce some of the relevant results required for our 
analysis. 


The 6 parameter in front of the the Chern-Simons term 
in Eq.(2.2) is taken to be 0 = ^ to ensure tha t th e 
statistics of the spins (hard-core bosons in Eq. (2.1)) 
are correctly transmuted to those of the fermions in Eq. 
(2.2). The Chern-Simons term on the kagome lattice can 
be written as 


=^5 + 45 

41 = / dt'^Ao{x,t)Ji{x - y)Ai{y,t) 

41 = -^ / dtJ2Mx,t)Kij{x - y)Aj{y,t) 


The first term 4^1 Eq.( |2.5[ ) is the flux attachment 

term that relates the density at a site on the lattice to the 
flux in its corresponding plaquette. Eor the case shown 
in Fig [2 the explicit expression for this term is given as 


J“(/c) = (l,-l,l,-e-*'=Ue-*'=Ul) 

J'(fc) = (0,e-*'=U-I,l,0,0) (2.6) 

J=(fc) = (-e-*'=U 0,0,0,-1,1) 


where ki and ^2 are the Fourier components along the ei 
and 62 directions of the unit cell shown in Figj^ These 
choices ensure that the fermion density n(x) (at a site x 
of the kagome lattice) is related to the gauge flux B{x) 
on the adjoining plaquette by the constraint equation 
n(x) = 6B{x) as an operator identity on the Hilbert 
space. 


The second term in Eq.(2.5) establishes the commuta¬ 
tion relations between the different gauge fields on the 
lattice and it is the structure of the Kij matrix that en¬ 


sures that the fluxes commute on neighboring sites. This 
condition is crucial to being able to enforce the flux at¬ 
tachment constraint consistently on each and every site 
of the lattice. The explicit expression for the Kij matrix 
is given as 



/ 

0 

-1 

1 

-S2 

Si -|- S 2 ^ 

+ S 2 ^ ^ 



1 

0 

i-sr' 

-S2 - Sj"^ 

Si 

-1 

1 


-1 

Si - 1 

0 

1 - S2 

Si 

-1 

2 



Si -|- S 2 ^ 

s^^ - 1 

0 

SlS^^ 




-S2 - Si ^ 



-S2SC^ 

0 

1-sr' 


1 

1 - S2 

1 

1 

-S2 

Si - 1 

0 / 


where Sj are shift operators along the two different di¬ 
rections (ei and 62 ) on the lattice i.e. Sjf{x) = f{x + ej) 
as shown in Fig[Tl 


III. THE AAZ MODEL WITH A CHIRALITY 
BREAKING FIELD 


Next, we will consider the effects of adding a chirality 
breaking term to the Heisenberg Hamiltonian in Eq.(2.1). 
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A system of spin-1/2 degrees of freedom on the kagome 
lattice with a chirality breaking term as its Hamiltonian 
was considered recently by Bauer and coworkersP^ Us¬ 
ing finite-size diagonalizations and DMRG calculations, 
combined with analytic arguments, these authors showed 
that the ground state of this system with an explicitly 
broken time-reversal invariance is a topological fluid in 
the universality class of the Laughlin state for bosons at 
level 2 (or, equivalently, filling fraction 1/2). Here we will 
examine this problem (including the XXZ Hamiltonian) 
and find that the ground state has indeed the same uni¬ 
versal features found by Bauer and coworkers, and by us 
in the 1 /3 plateau.!^ 

The resultant Hamiltonian is given as 

Htot = HxXZ + Hcl, — /lext ^ /g 


where Hxxz is the XXZ Heisenberg Hamiltonian in 
(2.1). The chirality breaking term is given by 


Hch = h'^Xijki^) = h'^Si- {Sj X Sk) 


A 


A 


(3.2) 


where Xijki^) is the chirality of the three spins on each 
of the triangular plaquettes of the kagome lattice and 
the sum runs over all the triangles of the kagome lattice. 
Recall the important fact that each unit cell of the of the 
kagome lattice contains two triangles. 

In order to use the flux attachment transformation, it 
is convenient to express the spin operators and 
in terms of the raising and lowering S'^ and S~. As an 
example, one can re-write the chirality term on a trian¬ 
gular plaquette associated with site b (shown in Fig[^ as 
follows 

X6 =Sa ■ {S, X Sb) 

=l\ -s-s+s^b+Sts- 

^ , , (3.3) 

I Q— QZ C+ C+ QZ Q— ^ ' 

-s^,s-s+ + s^,sts^^ 


where the subscripts a, b and c label the three corners of 
a triangular plaquette in Fig[2 

As shown in Ref. [5] (and summarized in Section 
the raising and lowering spin operators S^ are in¬ 
terpreted as the creation and destruction operators for 
bosons with hard cores, and S^ operators are simply 
related to the occupation number n of the bosons by 
S^ = I — n. Under the flux attachment transforma¬ 
tion, the hard core bosons are mapped onto a system of 
fermions coupled to Chern-Simons gauge fields (residing 
on the links of the kagome lattice). The boson occupa¬ 
tion number at a given site is mapped (as an operator 
identity) onto the gauge flux in the adjoining plaquette 
(in units of 27r). 

It is the straightforward to see that the chirality term 


gets mapped onto an additional hopping term on the links 
of the kagome lattice which carries a gauge as an extra 
phase factor on each link determined by the fermion den¬ 
sity on the opposite site of the triangle. As a re sult, only 
the fermionic hopping part of the action in Eq.(2.3) gets 
modified, and the interaction part and the Chern-Simons 
part are unaffected. 

Putting things together we get an effective fermionic 
hopping part that has the form 


{iDo + fi)tl;{x) 

- Y, J(a;(a)) + /i-c) 

{x,x') 

(3.4) 

where once again x and x' are nearest neighbor sites and 
X(^a) refers to third site on the triangle formed by sites x 
and x'. The subscript (a) refers to the sub-lattice label. 
The expressions of J and on each sub-lattice can be 
written as 


SF{tp,i’*,Af,) 


-iv 




1 

2 


'J 2 -h/l 2 



=tan 


h 

J 



(3.5) 


Hence, we have expressed the effects of the chirality term 
in terms of a modified hopping strength and an ad¬ 
ditional gauge field {cj)(^a){x)) on each of the links of the 
lattice. In the limit that ft, = 0, we just have the original 
gauge fields and in the other limit with J — 0 each link 
has an additional contribution of ((/(„)(a:) = ±|). 


IV. MEAN-FIELD THEORY 


In this section, we will set up the mean-field expres¬ 


sions for the fermionic action in Eq.(2.3l and Eq.(3.4|. 


The basic setup here is very similar to the situation de¬ 
scribed in our earlier work,l^ but it has been modified 
to account for the addition of the chirality term in this 
paper. 

Using the flux attachment constraint imposed by the 
Chern-Simons term (n(x) = 9B{x)), the interaction term 
in Eq.(2.3) can now be re-written as follows 


(A^) = AJ ^ Q - 0B{x)^ Q - 9B{x')^ ( 4 . 1 ) 


The interaction term has been expressed purely in terms 
of gauge fields. Hence, the resultant action after the flux 
attachment transformation is quadratic in the fermionic 
fields. Integrating out the fermionic degrees of freedom 
gives rise to the below effective action just in terms of 











6 


the gauge fields 

SesiAfj,) = -itT\n[iDo + fj, - i?hop(^)] , , 

+ 'S'int(^/i) + (^Scsi^fj.) 

where the hopping Hamiltonian H^hop(^) is (in matrix 
notation) 

^hop = ^ \x) {x'\ + h.c} (4 3) 

{x,x') 


where the above sum runs over all nearest neighbors x 
and x' . The gauge field Aj{x) refers to the hopping term 
required to go from point x to x' on the lattice. The 
term J(a) and are as defined in Eq.(3.5) with (a) 
once again referring to the sub-lattice index. In the above 
expression (a) would correspond to the third site in the 
triangle formed by nearest-neighbor sites x and x'. 


Now the mean-field equations can be obtained by ex- 
tremizing the action in Eq.(4.2) w.r.t. the gauge fields 


SA^ 



= 0 


(4.4) 


Differentiation with respect to the time components Aq 
yields the usual equation relating the density to the flux, 

{n{x)) = ^{B{x)) (4.5) 


which implies that the flux attachment is now enforced 
at the mean-field level. The average density can be ex¬ 
pressed in terms of the mean-field propagator by 


{n{x,t)) 


SSf 

6Ao{x,t) 


—iS{x,t',x,t) (4.6) 


where Sp refers to just the fermionic part of the action 
(i.e. the hopping part) and S{x,t]x',t') is the fermion 
propagator in an average background field A^(x,t)- 


Differentiation with respect to the spatial Ak compo¬ 
nents yields an expression for the local currents, 


, - , J l<Scs \ J JSl... 


(4.7) 


Here too, we can express the average current in terms of 
the fermionic action in the usual manner 


=(- 


SSf 


dAk{x,t) 


(4.8) 


We will look for uniform and time-independent solu¬ 
tions af these equations. Under these conditions the 



FIG. 2: Magnetic unit cell at half filling, a, b and c 
label the different sub-lattices in each of the unit cells. 
The gauge fields now have an additional label to 
indicate the unit-cell they belong to. 


mean-field equations for the currents, Eq.(4.7), becomes 


{jk{x)) =0d^°‘Aoaix) 

- 2JA02(-l)fc - fkS- - (1 - fk)B^] 

(4.9) 

with /fc = 1 when k = 1,5,6 and fk — 0 when k = 
2,3,4. In the above expression, we have also fixed the 
average fluxes on each sub-lattice (i.e. the fluxes on all 
sub-lattices of a particular type are the same), a is the 
sub-lattice index, and 



/ 

1 

0 

-S2 





-1 


0 



dka ^ 


1 

■52^ 

-1 

1 

0 

0 


(4.10) 


s 


0 

-1 




\ 

-1 

0 

1 



where si and S 2 are the 

same shift operators discussed 


after Eq.(2.7). 


A. Mean-fleld ansatz for XXZ model 

We will begin by studying the case with of the XXZ 
Heisenberg model i.e. we set to zero both the chirality 
coupling h and the external magnetic field, hext- This 
translates to the case of half-filling in the fermionic lan¬ 
guage. At 4 filling, the average density within each unit 
cell is given by 


^ {{na) + (n,) + (n,)) = ^ (4.11) 


where a, b and c refer to the three sublattices. This gives 
an average flux of tt in each unit cell which implies that 
the magnetic unit cell consists of two unit cells as shown 
in Fig[^ 

In the absence of the chirality term we will primarily 
look for mean-field phases that are uniform and time- 
indepe nden t, and h ave zero currents, i.e. {jkix,t)) = 0 
in Eq.(4.5) and Eq.(4.7). The flux attachment condition 
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can be imposed as follows on each of the sub-lattices 
{na{x)) =0{Ba{x)) = i - Ai - A 2 
{nb{x)) =0{Bb{x)) = ^ -f Ai (4.12) 

{nc{x)) =0{Bc{x)) = ^+A 2 


where Ai and A 2 are two parameters that will be chosen 
to satisfy the mean-field self-consistency equations. The 
fluxes in Eq.( 4.12| can be achieved by the below choice 
of gauge fields in Fig[^ 


A\=0 Al=0 

Al =pi Al= Pi 

yfi = 0 Al = 0 

^1 = 0 Al = 0 

Al = -P2 Al = -P2 + Stt 

^ 1=0 AI = 3tt 


(4.13) 


where pi = tt + 27rAi and p 2 = tt + 27rA2. 

With these expressions for the densities, the mean-field 
equation (Eq.( |4.9[ )) can be satisfied by the below choices 
for the temporal gauge fields 


^o,a — 2A(Ai -|- A 2 ) ^0,6 — —2AAi Aq^c — — 2 AA 2 

(4.14) 

Using this mean-field field setup, we find two regimes at 
the mean-field level. 


1. XY regime 

In the XY regime, j < 1, we find that Ai = A 2 = 0 
is the only solution that satisfies the self-consistency 
condition. This leads to a state with a flux of tt in 
each of the plaquettes. We will represent this as the 
((i?a), {Bb), {Be)) = ( 7 r, 7 r, 7 r) flux state. This state has a 
total of six bands, shown in Fig. (the top two bands 
are double degenerate). At half-filling the bottom three 
bands are filled giving rise to two Dirac points in the spec¬ 
trum, crossed by the dotted line in Figj^ which indicates 
the Fermi level. See SeclV]for details. 

At the mean-field level this spectrum is equivalent to 
the gapless U{1) Dirac spin liquid state that has been 
discussed in previous works.^^^^^ We notice, however, 
that there are other works that favor symmetry breaking 
states but with a doubled unit cell and a flux of tt in 
each of the plaquette^. The state we find could survive 
when fluctuations are considered giving rise to one of the 
above states. Alternatively, fluctuations could also open 
up a gap in the spectrum leading to an entirely different 
phase. In this paper, we will only analyze the gapless 
states at a mean-field level. 


En 



FIG. 3: (Color online) Mean-field spectrum in the XY 
regime at half-filling, showing the two Dirac points. The 
dashed line indicates the Fermi level. The top band is 
doubly degenerate. These plots are made along the 
ky = —^ line in the Brillouin zone (along which the 
two Dirac points lie). 


2. Ising Regime 


For 7 ^ Ij non-vanishing values of Ai and A 2 are 
required to satisfy the mean-field consistency equations. 
The solution with the lowest energy has the form Ai = 
—A 2 ^ 0. This solution shifts the mean-field state away 
from the pattern (tt, tt, tt) for the flux state, and opens 
up a gap. 

The Chern number of each of the resulting bands can 
be computed by using the standard expressiorP^ in terms 
of the flux (through the Brillouin zone) of the Berry con¬ 
nections 


C=^ [ d^kX,y{k) (4.15) 

JBZ 

where Xij = diAj — djAj is the flux of the Berry connec¬ 
tion Ai = —i ('01 dki 10). Here |0) refers to the normal¬ 
ized eigenvector of the corresponding band. 

In the Ising regime the Chern numbers of the bands 
are 


Cl = 0 C 2 = 0 C 3 = 0 

C4 = 0 C5 = 0 Cg = 0 


(4.16) 


This implies that in the Ising regime, the total Chern 
number for the filled bands is 0. This means that we are 
left with the original Chern-Simons term from the flux 
attachment transformation. In this regime, the fermions 
are essentially transmuted back to the original hard¬ 
core bosons (spins) that we began with and our analysis 
doesn’t pick out any specific state. 














FIG. 4: Currents induced by the chirality term in each 
of the triangular plaquettes. The currents are indicated 

by j- 


B. Mean-fleld theory with a non-vanishine chirality 
field, h^O 


In this section, we will turn on the chirality term. 
Looking at the doubled unit cell in Fig. there are 
four corresponding chirality terms (within each magnetic 
unit cell) which can be written as 

Xtjkix) = Xbi{x) + Xb 2 {.x) + Xci{x) + Xc2(x) (4.17) 


Now, we haye to account for the additional contribu¬ 
tions from h in Eq.(3.5|. Importantly, the added con¬ 
tribution to t he g auge fields due to a non-zero yalue 
4'{a)ix) in Eq.(3.5) will giye rise to additional fluxes and 
shift the state away from the (tt, tt, tt) flux state obseryed 
in the XY regime section in Sec |IV A 1| Notice that, 
if we were stay in the (tt, tt, tt) flux state, this would 
imply that the ayerage density (n(a)(x)) = ^ at eyery 
site. In this situation the expectation value of the chiral¬ 
ity operator automatically vanishes due to the relation 
^^{a)) ~ \ ~ this situation the chirality 

term would never pick up an expectation value at the 
mean-field level and time reversal symmetry would re¬ 
main unbroken. Hence, in a state with broken time re¬ 
versal invariance the site densities cannot all be exactly 
equal to 

The fluxes in each of the plaquettes also gets modified 
due to the contribution from (f^^x). The effective flux at 
each of the sublattice sites is now given as 


clockwise manner. Hence, we will choose an ansatz on 
each of the different links as seen in Fig. [4| T he mean- 
field equations for the current terms in Eq.(4.9) can now 
be satisfied by the below choice of gauge fields 


^ 0 ,a —2AJ(Ai -|- A 2 ) + j 
Ho,fa=-2AJAi-J (4.19) 

Aq^c = — 2A J A 2 — j 


In the above equations the effect of the chirality term di¬ 
rectly enters in the form of a current. Now, we will pro¬ 
ceed to look for mean-field phases that self-consistently 
satisfy the mean-field equations in Eq.(4.5l and Eq.(4.9) 
as well as constraints set by Eq.(3.5) and Eq.(4.I8). Once 
again, we will analyze the cases of the XY and Ising 
regimes separately. 


1. The XY regime 

In the XY regime, j < 1, we had the (tTjTt, tt) flux 
state which was gapless and had two Dirac points (see 
Fig. §. Here, we find that even for small values of h, 
there exist solutions with Ai = A 2 7 ^ 0. This shifts the 
state away from the (tt, tt, tt) flux state and opens up a 
gap in the spectrum as shown in Fig. 

The values of the mean-field parameters for a few dif¬ 
ferent values of the field h (the strength of the chirality 
breaking term) are shown in Table |lj A plot of the mean- 
field spectrum for the specific case of /i = 0.05 J is shown 
in Fig. As the value of h is increased from 0, the av¬ 
erage flux on each of the triangular plaquettes decreases 
from TT —The corresponding flux in the hexagonal 
plaquettes goes from tt —>■ 27r. In the limit of a strong 
chirality term, one would expect to get a state with flux 
of 27r in each hexagonal plaquette, and a flux of ^ in each 
of the triangular plaquettes. We will refer to this as the 
(27r, 1 ) flux phase. The values of the energy gap and 

the expectation values of the chirality operator are also 
shown for the different values of h in Table |l] (The en¬ 
ergy gaps essentially measure the gaps between the Dirac 
points.) 

The Chern numbers of the six bands with the chirality 
breaking field h turned on are 


{Ba) =7r - 27rAi - 27rA2 -\-2{(j)a + (j)b + (t'c) 

{Bb) =TT + 2-kXi- 4>b +4>c) (4-18) 

{Be) =7r -1- 27rA2 - {(j)a + 4>b + 4>c) 


The above fluxes still ensure that we in the half-filled 
case. 

In order to accommodate such a flux state, we also have 
to allow for non-zero currents in the mean-field state in 
Eq. (4.9). As a result we will consider an ansatz with 
{jk{x, t)) 7 ^ 0. The chirality terms in the Hamiltonian 
go across each of the triangular plaquettes in a counter- 


Ci — -fl C2 — —I Gs — -fl 
Ci = -fl G5 = —I Gg = —1 


(4.20) 


As we are still at half-filling, the bottom three bands 
must be filled, leading to a total Chern number of the 
occupied bands of Gi + G 2 + G 3 = -1-1. This along 
with the original Chern-Simons term from the flux at¬ 
tachment transformation is expected to give rise to an 
effective Chern-Simons term with an effective parameter 

^eff = Of + Ocs = X—I- X— (4-21) 
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h 

J 

II 

Ai — A 2 

Eg 

(x) 

0 

0.500 

0 

0 

0 

0.05 

0.460 

-0.040 

0.2286J 

0.000782 

0.1 

0.385 

-0.115 

0.5518J 

0.001064 

0.5 

0.300 

-0.200 

0.7351J 

0.002149 

1 

0.275 

-0.225 

0.7638J 

0.002642 


TABLE I: Approximate values for the mean-field parameters for different values of chirality (h) for j = 1- Here Eq 
denotes the energy gap in units of J, and (y) is the expectation value of the chirality operator. As the chirality term 
gets stronger, the average density on each of the triangular plaquettes approaches 0.25, and the density on the 

hexagonal plaquettes approaches 1 . 


En 



FIG. 5: (Color online) The chirality term opens up a 
gap in the spectrum (see Fig[^ The above plot is for 
h = 0.05J. The plot is made along the in 

the Brillouin zone. 


A more detailed and rigorous computation of the above 
statement will be presented in a later section. Sec. |Vj 
where we will include the effects of fluctuations and show 
that the resultant continuum action is indeed a Chern- 
Simons theory with the above effective parameter. This 
result shows that in the presence of the chirality term, we 
do obtain a chiral spin liquid. Such a state is equivalent to 
a Laughlin fractional quantum Hall state for bosons with 
a spin Hall conductivity The state obtained 

here has the same topological properties as the state that 
we founcP in the magnetization plateau at m = ^. 


In the limit that we only have the chirality term, we 
have that J = 0. Now = ±f in Eq.(3.5|. We can 
now look for solutions such that 



4^b — — 


TT 

2 


(4.22) 


Here, we again recover the (27r, f, f) flux state as ex¬ 
pected. The resulting chiral state is the same found by 
Bauer and coworkers in Ref. |25j . 


2. Ising regime 

In the Ising regime, j > 1, the Heisenberg model gave 
rise to a state that was gapped and a vanishing Chern 
number, as shown in Sec. |IV A 2| Here a small chirality 
term would not affect the mean-field state as long as it is 
weak enough. In order to see the chiral spin liquid state 
obtained in the XY regime in Sec. |IVB 1[ one would 
need a strong enough chirality term to close the Ising 
anisotropy gaps and to open a chiral gap so as to give rises 
to states with non-trivial Chern numbers. Hence, the 
state here would be determined based on the competition 
between the anisotropy parameter A and the strength of 
the chirality parameter h. 


C. Combined effects of a chirality symmetry 
breaking term and an external magnetic field 


So far, we have primarily focused on the case of half¬ 
filling and hence in the absence of an external magnetic 
field, hext = 0- Now we will briefly consider the scenario 
when the external magnetic field is present, /text ^ 0, 
in Eq. (3.11 or, equivalently, that we are at fermionic 
fillings other than ^ in the XY limit. This will allow us 
to connect our recent results on a chiral spin liquid phase 
in a magnetization plateau with the chiral state arising in 
the presence of a chirality symmetry breaking held. The 
mean held theory we discuss here has points of contact, 
including the role of Chern numbers, with a classic paper 
by Haldane and Arovas.l^ 

In the previous section we noted that the main effect 
of adding the chirality symmetry breaking term to the 
mean-held state was to shift the huxes on each of the 
sub-lattices. We began with a (tt, tt, tt) hux phase for 
the Heisenberg model and it was modihed to a (27r, f, f) 
hux phase in the presence of a strong chirality term. Es¬ 
sentially the chirality term shifted the huxes from the 
triangles to the hexagons. Using this analogy, we will 
now look for similar hux phases at other hllings. In the 
presence of a strong chirality term, we will consider hux 
phases where the hux in maximized in the hexagons and 
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minimized in the triangles at different fillings. 

In the absence of the chirality term, we have a uniform 
flux phase with {Ba) = (Bb) = (bJ) = (j) = 27r| with 
p,q ^ h. When, we turn on the chirality term, we expect 
the fluxes from the triangles to shift to the hexagons. 
Hence, we have 


{Ba) =(/) + 25 

{Bb)=^-6 (4.23) 

(Be) =4>- 5 

so that the total flux in each unit-cell is still the same. 
Such a flux state can be realized by the below choice of 
gauge fields 

Hi(5f)=0, A2{x) = (j) — 5, 

A:i{x) = 0, AAS) = 0, (4.24) 

Ab{x) =—(p + 6 + Aq{x)=3(Pxi 

with X = {xi,X 2 )- xi and X 2 are the coordinates along 
the Cl and 62 directions in Fig. [^respectively. The fluxes 
on each plaquette range between 0 and 27r, which trans¬ 
lates to having a site filling between 0 and 1. Hence, we 
set 5 = min(^, TT — ^). Using this choice, one can plot 
the Hofstadter spectrum in the limit of a strong chiral¬ 
ity term. In Fig. [^we plot the Hofstadter spectrum for 
the case with j = 5. The bottom solid line indicates the 
Fermi level (all the occupied states) and the top solid line 
indicates the next excited state available. 

At most fillings the total Chern numbers of all the oc¬ 
cupied bands is —1. This would lead to a Chern-Simons 
term with pre-factor — ^ and such a term would be ex¬ 
pected to cancel when combined with the original Chern- 
Simons term from the flux attachment transformation, 
which also has a pre-factor ^. The exceptions are at the 
fillings (n) = g) represented by vertical jumps in 

the solid lines in Fig. [^ At these fillings, the total Chern 
number of all the filled bands is different from —1 and 
lead to an effective Chern-Simons term. The resulting 
magnetization plateaus and their corresponding Chern- 
numbers are summarized in Table El 


(n) 

Chern No. 

m 

1 

6 

-hi 

2 

3 

1 

3 

+1 

1 

3 

4 

9 

+2 

1 

9 

1 

2 

-hi 

0 


TABLE H: Magnetization plateaus obtained in Fig. [^ 
and their Chern numbers. At these fillings the 
Chern-Simons terms do not cancel out and the system 
is in a chiral spin liquid. 


have also been previously obtained in the absence of the 
chirality symmetry breaking term, h = 0 and /igxt A 0 
in Ref. [5]. It is apparent that these plateaus survive 
in the presence of the chirality symmetry breaking term. 
Additionally we observe two other plateaus at fillings | 
and i. The plateau at ^ filling is the same one that was 
observed in the previous sections for the case with no 
magnetic field (see Sec. IVBlI. 


The plateau at | filling has a magnetic unit cell with 
three basic unit cells. This gives rise to a total of nine 
bands of which four are filled. The Chern numbers of 
each of the nine bands in the mean-field state are 


C\ — —1 C 2 — -\-2 C 3 — -\-2 

Ci = -1 C5 = -1 Cg = +2 (4.25) 

C'j = —1 Cg = —1 Cg = —1 


The Chern numbers of the four filled bands (Ci, C 2 , C 3 , 
and C 4 ) add up to +2. Again this result will combine 
with the Chern-Simons term from the flux attachment 
transformation leading to an effective Chern-Simons with 
an effective spin Hall conductance of ai,, = I. This result 
is also summarized in Table [H| In Ref. [5] we identified 
this state as having the same topological properties as 
the first state in the Jain sequence of fractional quantum 
Hall states of bosons. 



FIG. 6 : (color online) Hofstadter spectrum for h = 5J. 
The cc-axis plots the mean-field fermion density (n) and 
the y-axis plots the XY energies. The bottom solid line 
indicates the fermi level and the top solid line 
corresponds to the next excited state available. The 
vertical jumps in the figure correspond to possible 
magnetization plateaus. 


The magnetization plateaus at filling fractions ^ and | 
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D. Chiral Spin Liquids with Dzyaloshinski-Moriya 
Interactions 


In this Section, we consider the effects of a 
Dzyaloshinski-Moriya term (instead of the chirality term) 
on the nearest neighbor XXZ Heisenberg Hamiltonian in 
Eq. (2.1 1 . The Dzyaloshinski-Moriya term is written as 


Hdm = Jdm ^ z • {Si X Sj) 


(4.26) 


where the sum runs over nearest neighbors in each 
triangle in a clockwise manner. As an example the 
Dzyaloshinski-Moriya term in a triangle associated with 
site bi in Fig. [^can be written as 

Hdmm = ("5'+(x)S'j7^(f) - S'+(x)5'^Jf) 


+S+^{x)S^^{x) - S+^{x)S^^{x) 
{x)S-^{x) - {x)S-^{x) 


(4.27) 

Clearly this term breaks time reversal so we expect that 
we may be able to find chiral phases. 


From the form of Eq.(4.27), we can now readily apply 


the flux attachment transformation just like we had for 
the case of the chirality term. As a result the parameters 
in Eq.(3.5l now get modified as 






(l>(a){x) =tan ^ j Q - n(a){x)^ + 


Jdm\ 

~ J 


(4.28) 


In this section, we will set the chirality symmetry break¬ 
ing term to zero, h = 0 . 

Two separate regimes have to be considered. 


1- Jdm < 1.7J 


Recall that in the XV regime the Heisenberg model 
gave rise to the (7r,7r,7r) flux state which is gapless and 
has Dirac points (Fig. [^. Treating the Dzyaloshinski- 
Moriya term as a perturbation, we find that this term also 
opens up a gap in the (tt, tt, tt) flux state as can be seen in 
Fig.0 But a the resultant state obtained still has a flux 
of TT in each of the plaquettes. For the situation shown 
in Fig[^ the energy gap is 0.I366J. This is an important 
difference between the effects of adding the chiral term 
and the Dzayloshinkii-Moriya terms, since the chirality 
term shifted the fluxes on each plaquette away from tt 
whereas the Dzyaloshinski-Moriya term does not. 

The Chern numbers of the six bands in the presence of 



FIG. 7: (color online) Spectrum with Jbm = 0.05 J. 
The Dzyaloshinski-Moriya term opens up a gap in the 
spectrum of the Heisenberg Hamiltonian (see Fig. |^. 


a small Jdm term are 


Cl — -l-I C2 — —1 C3 — -|-I 
C 4 = +l C 5 = -1 Ce = -1 


(4.29) 


Once again, we find that the total Chern number of all 
the filled bands is Ci -I- (72 -|- Cs = -l-I. This would again 
lead to a fractional quantum Hall type phase with a^y = 
i, just as we had observed in the case with the chirality 
term. 


E. Jdm ^ 1.7J 


For larger values of the Dzyaloshinski-Moriya param¬ 
eter, namely for TaM. > x.7, the Chern numbers of the 
bands again get rearranged and the chiral phase no longer 
survives, as shown below 


Cl — —I C 2 — +1 C 3 — —I 
C 4 = —I C 3 = -|-1 Cq = -|-I 


(4.30) 


In the limit that we only have the Dzyaloshinski- 
Moriya term, i.e. J = 0, the values of all = | in 
Eq.(4.28). In this case the values of the mean-field pa¬ 
rameters that satisfy the consistency equations, Eq.(4.5) 


and Eq. (|4.9|), are Ai = A 2 = — 4. Hence, in the pres¬ 


ence of only t' 


le Dzayloshinski-Moriya term, we again end 


up in the (tt, tt, tt) flux state that was observed in the XY 


regime of the Heisenberg model in Sec IV A 1 


F. Dzyaloshinski-Moriya term with an nniform 
magnetic field, /lext 7 ^ 0 


Finally, we will also consider the effects of the 
Dzyaloshinski-Moriya term in the presence of an uniform 
external magnetic field /lext 7 ^ Oj just as we had done 
for the chirality terms in Sec. |IV C We will once again 
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(n) 

Range of values (in J) 

Chern No. 

m 

1 

6 

0 < Jdm ^ 0.35 

+1 

2 

3 

2 

9 

0 < Jdm ^ 0.3 

+2 

5 

9 

1 

3 

0 < Jdm ^ 0.8 

+1 

1 

3 

4 

9 

0.05 < Jdm < 0.6 

+2 

1 

9 

1 

2 

0 < Jdm ^ 1-7 

+1 

0 


TABLE III: Approximate values of Jdm for which we 
observe the plateaus at the mean-field level. This table 
also lists the corresponding Chern numbers and their 
magnetizations, m. 


field term again. 


FIG. 8: (Color online) Hosftadter spectrum in the XY 
limit for the case Jdm = 0.3J. (The x-axis plots the 
mean-field fermion density (n) and the y-axis plots the 
energies of the XY model.) The bottom solid line 
represents the Fermi level (all the filled bands) and the 
top solid line represents the next excited energy state 
available. The plateaus correspond to jumps in the solid 
line. We see two additional plateaus at densities 
(n) = |, I for certain values of Jdm- 


focus on the XY limit where the mean-field equations 
are simpler due to the absence of the interaction term, 
i.e. A = 0. We will look for states that are uniform, 
time-independent and don’t have any currents. 

This scenario is very similar to the case of the inte¬ 
ger quantum Hall effect with non-interacting fermions in 
the presence of a (statistical) gauge field. This approach 
was also used by Misguich et. al. in their studies on the 
triangular lattice.!^ More recently, we carried out a sim¬ 
ilar analysis on the kagome lattice with an XY nearest 
neighbor Heisenberg model.^ Here, we will perform the 
same analysis, but with the Dzyaloshinski-Moriya term 
added to the XY nearest neighbor Heisenberg model. 

Once again, we find a few different plateaus as can 
be seen in the Hofstadter spectrum in Fig for Jdm = 
0.3J. The vertical lines in the figure correspond to the 
magnetization plateaus. The range of Jdm values for 
which we observe the above plateaus is shown in Table 
mil The table also lists the total Chern numbers of all 
the filled bands at each of the plateaus as well as the 
corresponding magnetization. 

This concludes our mean-field analysis into the various 
possible magnetization plateaus. We will now proceed to 
consider the effects of fluctuations on the mean-field state 
when a small chirality term was added to the Heisenberg 
model in the XY limit. This was the situation discussed 
in Sec |IVB[ For the rest of the paper, we will not consider 
the Dzyaloshinski-Moriya term or the external magnetic 


V. EFFECTIVE FIELD THEORY 


In this section we return to the case of the nearest 
neighbor Heisenberg model in the presence of a small chi¬ 
rality term. In Sec. |IVB[ it was shown that the addition 
of the chirality term opened up a gap in the mean-field 
spectrum and lead to a state with non-trivial Chern num¬ 
ber. We will now expand the fermionic action around 
this mean-field state and consider its continuum limit. 
This process will allow us to go beyond the mean-field 
level and consider the fluctuation effects of the statistical 
gauge fields. The analysis presented here is analogous to 
the one presented in our earlier work.l^ As a result we will 
only write down the relevant expressions for the current 
scenario. 


In Sec. |IVB I[ we found that in the absence of the 
chirality term the spectrum was gapless with two Dirac 
points and that the addition of the chirality term opened 
up a gap at the Dirac points. These two Dirac points in 
the mean-field phase were located at the momenta K = 


± 


. The fermionic degrees of freedom on each 
site can be expanded around each of the two Dirac points 
using the following expansions on each sub-lattice 


/V’alA 
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(5.1) 
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where in r refers to the Dirac species index and the 
label a refers to the spinor index within each species. 
tpai,iJbi,il’ci,tpa 2 ,tpb 2 , 4 ’c 2 refer to the original fermionic 
fields on the different sub-lattices sites in the mean-field 
state at half-filling as shown in Fig[^ 

Now we will include the fluctuating components i.e. 
we will expand the statistical gauge fields as follows 
= {A^) -I- 5A^. The mean-field values of {A^) are 
the same as those given in Sec. |IV B| From now on, we 
will primarily focus on the fluctuating components. In 
order to simplify the notation, we will drop the 5 label 
in the fluctuating components i.e. all the gauge fields 
presented beyond this point are purely the fluctuating 
components. 


A. Spatial fluctuating components 

First, we will begin by looking at just the spatial fluc¬ 
tuating components. Furthermore, we will also expand 
all the spatial fluctuating components in the magnetic 
unit cell in Fig. [^in terms of slow and fast components. 
This will allow us to treat the slow components as the 
more relevant fields. 

The fields along the ei direction (in Fig. can be 
expanded as 

+ -Af) 

+ + Af) 

aI^^(a^-a{^-a{^ + a{^) 

Similarly, the fields along the e 2 direction can be written 
as 


A\^^ 
^ 2 

Al='^ 
® 2 

^ 2 

® 2 


-i/l 4-^ A _ _ 4^2 _L 4/3 

-h ^ A±y Ai.2 ^2 ' ^2 

-lA. + ^A, + Ai^+Ai^-Ai^ 

-I A 4-—A _ 4ti + 4t2 _ / 1/3 


(5.3) 


In the above expressions A^^'^ refer to the fluctuating 
components along the different links of the unit cell in 
the mean-field state (1 and 2 refer to the two unit cells 
in the magnetic unit cel shown in Fig. [^. The slow 
components are represented by A^ and Ay and A^, ^ 

A{'^ and A{^ are the fast fields along the different spatial 
directions. 


B. Temporal fluctuating components 


Similarly, the fluctuating time components can also be 
expanded in terms of slow and fast fields as follows 


Ao,.i =ao (Ao + A^°i + A^“i + 3V2{AI^ - A^^)) 

Ao,bi =0-0 ^Ao — Aq^° — 3Aq^ j 

Ao.ci =ao (^Ao - AI°^ -t 3\/2-V3AI^ - 3\J2 + VsA^^ 

Ao,a 2 =ao (Ao + A^ + Ar - 3V2iA{,^ - A^^)) 

Ao,i,2 =oo ^Ao — Aq^° -|- 3Aq®^ 


Ao.ci =00 (^Ao - A^°^ - 3^/2-V 3A^^ -t 8^/2 -t ^Af ^ 

(5.5) 

where al, &1, cl, a2, b2 and c2 again refer to the dif¬ 
ferent sub-lattice indices in the mean-field state in Fig[^ 
The only temporal slow component is Aq. All the other 
fields with super-script / refer to the fast fields. The pre¬ 
factors and constants in Eq.(5.5) are chosen to make the 
notation and computation below easier. 


Using Eq.(5.1), the mean-field action with the choice of 
the mean-field gauge fields in Sec IV B in the continuum 
limit becomes 


Sf,mf 


— ‘rn)'^r 


(5.6) 


where A' = ^'*70 with Af = (5'J, A'f, A'f, We are 

using the slash notation $ = = 70^0 — lidi with the 

Minkowski metric g^i,. The gamma matrices act on the 
upper or spinor index (a) in fkj* and are given by 


Finally, the fields along ci -I- 62 directions can be ex¬ 
pressed as 

4 =? + ^Ay + ^Af - Ai^ - At + Af) 

Ag = Y ( 2 "^^ ~ 

^ ^ (5A) 

Al =1 ( 1 a. + fAy + |a^ - Af - At + At) 

Al (^A. + ^Ay + |a^ - A^ + At - At) 


7o = 0-3 7i = *cr 2 72 = * 0-1 (5.7) 


Importantly the mass terms m are the same for both the 
Dirac points and are given as 


m = lim 

ao->-0 


-J^(3 + y3) 

9ao 


0.5258 

ao 


A > 0 (5.8) 


Hence, the masses m are positive for both Dirac points 
as the value of A < 0 from the mean-field analysis (as 
shown in Table |T| . 

The resulting action for the spatial fluctuating compo- 
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nents becomes 

^^siow ^ J ^3^ + AyW724'} 

^ J d3a;|^A'^4'7i4'+ '^A/1'724' 

_ (Tlf _ ^/3 ^ ^/2 _ ^/3 _ ^/2 ^ 

(5.9) 

where we have absorbed some of the constant factors into 
the definitions of the fast fields to make the notation more 
convenient and the definitions of the gamma matrices are 
the same as in Eq. (5.7). 


action for the slow components becomes 
Sf, s\ow = / d^x'i - to] 4 ' 


(5.11) 


The resulting continuum action for the slow and fast 
fields become 

= J d^xl^Al^ (^7ir34-) + ylf (^72T34') 

+ aI^ (^70^3^) 


where is the covariant derivative. The 

fast components can be expressed as 

SFMt = J + 

+ + ^^^1'72T^4' (5.12) 

where 

0 . ( a {^ _ a (3 + 43 _ 43 _ 43 ^ 43^ 

(5.13) 


Eq.(5.11) and Eq. (5.12) can also be expressed in mo¬ 


mentum space as 


Sf = 




(2f 


(5.14) 


(5.10) 

where is the regular as Pauli matrix but acting on 
the species index r in 41^. Combining equations Eq.(5.6), 


Eq.(5.9) and Eq. (5.10), the total continuum fermionic 


with M = M + SM. 

The mean-field part M is given as 

y U p — m 

and the fluctuation part 5M is given as 


(5.15) 


y 0 -4-4>i+(^Af+ -A^^'jjs-A^^joJ 

(5.16) 


where (j)i is given in Eq.(5.13). 


The action in Eq.(5.14) is quadratic in fermionic fields 


and fermions can be integrated out to give an effective 
action in terms of just the fluctuating gauge fields. The 
resulting effective action becomes 


r 


ating components gives 


=lTr {m ^SMM ^SM 


tr [S{p)dM{q)S{p + q)6M{-q)) 


(5.19) 


S'eff = —jTrlnM 


(5.17) 


where M is defined in Eq.(5.16). Now we can expand M 


in terms of the mean-field part and the fluctuating parts 
as shown in Eq.(5.15) and Eq.(5.16). 


where the lower-cased ‘tr’ is a matrix trace, and S{p) = 
M{p)~^ is the continuum mean-field propagator pre¬ 
sented in Eq. (5.15), and it is given by 


S{p) = 


1 


p"^ — m 


f + m 0 
2 f 0 f + m 


(5.20) 


S^s = -fTr{ln {M + 5M)] 


(5.18) 


= -fTr{lnM}-fTr|ln(|l-tM | 

Expanding this action up to second order in the fluctu 


In the expansion of Eq.(5.18) we will only keep the most 


relevant (mass) terms (without derivatives) for the fast 
components. 

Similarly, one can also express the lattice version of 
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the Chen-Simons term and the interaction terms using 
the slow and fast fluctuating components. Combining all 
of the above terms, one can obtain the final continuum 
action. All the massive fields can safely be integrated out. 
This leaves us with just the Chern-Simons and Maxwell 
terms. The computation of this Feynam diagram is stan- 
dard and it is done in many places in the literature.^^^^^ 
To lowest order, after integrating out all the massive 
fields, the most relevant term is the effective Chern- 
Simons term , since it has the smallest number of 
derivatives, and is given by 

= {1 + y)J (5-21) 


where 0 = -^ from the original flux attachment transfor¬ 
mation and is the obtained from integrating out the 
fermions and is given as 


Of = ^(sgn(m) -F sgn(TO)) = ^ (5.22) 


as sgn(m) = -1-1 (m > 0 as shown in Eq.(5.8)). 


Hence, the Chern-Simons terms add up, and we get a 
state with spin Hall conductivity a^y = This state 
is equivalent to a bosonic Laughlin fractional quantum 
Hall state. This agrees and verifies our expected result 
obtained in Sec. IIVBI 

The Maxwell terms can be conveniently expressed in 
terms of the electric E and magnetic B fields as follows 


Sem 




(5.23) 


^ X = (24V3ao - . 

The computation in this section confirms our expec¬ 
tation and analysis used to determine the nature of the 
chiral spin liquid states using the mean-field theory ap¬ 
proaches in Sec. EYi 


VI. SPONTANEOUS BREAKING OF TIME 
REVERSAL INVARIANCE 


In the cases discussed so far in this paper, we be¬ 
gan with a (tt, tt, tt) flux state which , at the level of 
the mean field theory, has massless Dirac fermions, and 
showed that breaking the time-reversal symmetry ex¬ 
plicitly, by adding either a chirality term (Sec. IV B) 


or a Dzyaloshinski-Moriya term (Sec. IVD), led to a 


gapped state. We the showed, that quantum correc¬ 
tions led directly to a chiral spin liquid with broken time- 
reversal symmetry for arbitrarily small values of the chi¬ 
ral held h or the Dzyaloshinski-Moriya interaction Jdm- 
The existence of an explicit gap in the spectrum of the 
fermions was essential to this analysis. Furthermore, af¬ 
ter the leading quantum corrections are taken into ac¬ 
count, we found that the naive Dirac fermions of the 


mean-field theory became anyons (semions in the cases 
that were discussed in detail). This line of reasoning 
parallels the theory of the fractional quantum Hall effect 
where, at the mean held level, one begins with composite 
fermions fulling up effective Landau levels ,1^ whi ch tur n 
into anyons by virtue of the quantum corrections 

We now turn to the question of whether it is possi¬ 
ble to obtain a chiral spin liquid by spontaneous time- 
reversal symmetry breaking. This concept was formu¬ 
lated originally by Wen, Wilczek and Zee^ in the con¬ 
text of the J 1 — J 2 Heisenberg model on the square lattice, 
where a chiral ly-in variant Z 2 spin liquid appears to be fa¬ 
vored instead.l^I^ In this section we will show that ring- 
exchange processes on the bow ties (i.e. two triangles 
sharing the same spin) of the kagome lattice may favor 
the spontaneous formation of the chiral spin liquid if the 
associated coupling constant is large enough. Unfortu¬ 
nately, the critical value of this coupling constant that 
we obtain is much too large for the mean field theory to 
be reliable and, hence, we cannot exclude the possibility 
that other states may arise at weaker coupling. Nev¬ 
ertheless, it is an instructive excercise that shows that 
ring-exchange processes, if large enough, may trigger a 
chiral spin liquid on their own. 

In this section we explore of the possibility of break¬ 
ing this symmetry spontaneously. Numerical works 
have studied examples where such scenarios arise in the 
Heisenberg model on the kagome lattice in the presence of 
second and thi rd ne xt nearest neighbor Heisenberg terms 
or Ising terms,l^^^^ where they find suggestive evidence 
of a chiral spin liquid in certain regimes. Unfortunately, 
the flux attachment transformation summarized in Sec. 
[IT] cannot be applied to next nearest neighbor Heisenberg 
terms. However, we have examined the case in the pres¬ 
ence of just the next nearest neighbor Ising terms, using 
flux attachment methods and we do not find the chiral 
spin liquid observed in the numerical work.E^ 

As a result we consider the effect of adding a chiral 
term on a bowtie in the kagome lattice which is written 
explicitly as follows 


H^ = gY.[S.- {S, X Sk)] [S, ■ {Si X S,„)] 

{X} 

where the sum runs over all the bowties of the kagome 
lattice, i, j, and k refer to the indices of the up triangle 
and i, I and m refer to the indices of the down triangle, 
with i being the common site in the bowtie. 

The total Hamiltonian used in this section, can then 
be written as 


Htot —Hxxz + 

^Hxxz + ^ gX/^Xv ( 6 - 2 ) 

(A,V> 

where the Hxxz is the Hamiltonian of the Heisenberg 
antiferromagnet on the kagome lattice, with anisotropy 
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coupling A, defined in Eq. (2.1), and where XA and Xv 
are the chiralities over the up and down triangles (i.e. 
the sites of the two sublattices of the honeycomb lattice) 
and the sum runs over nearest-neighbor triangles of the 
kagome lattice (which correspond to the bowties). In 
what follows we will assume that we are either at the 
isotropic point or in the regime of XY anisotropy (easy 
plane), i.e. A < 1. 


We now note that the bowtie terms of the Hamiltonian 
in Eq.(6.1|, when expanded, can be expressed in terms 
of a ring-exchange term on the bowtie as follows 


f ^ i Sz ■ S, [{Sj ■ Si){Sk ■ Sm) - {Sj ■ Sn^){Sk ■ Si)] 

A ^ 

+ Si ■ Si [(Sj • Sm){Si ■ Sk) — {Si ■ Sj){Sk ■ Sm)] 

+S, ■ Sm [(S, ■ Sj){Sk ■ Si) - {Sj ■ S,)(S. • Sk)] I 

(6.3) 

Ring exchange terms have been known to give rise to ex¬ 
otic dimer states in Heisenberg antiferromagnets.^^Here, 
we will explore the possibility of such a term giving rise 
to a chiral spin liquid state. 


following effective action 
S = Sxxz 

+ J dt-^'^h{r,t)K-^{r,r')h{r',t) 

r,r' 

-b Seff[h(r,t),H^(r,t)] (6.5) 

where Seff [h(r, <), H^(r, t)] is the effective action of the 
fermions in a background chirality field h{r, t) (and which 
includes the lattice Chern-Simons term, as before). 

We can now carry out a mean-field approximation by 
extremizing the action of Eq.( |6.5[ ) with respect to the 
chirality field h{r,t), and to the gauge field Since 
we are working at zero external magnetic field, the mean 
field state for the gauge held is just the (tt, tt, tt) hux state 
and, hence, in the absence of any other interactions, we 
will naively have two species of massless Dirac fermions 
(as discussed in Sec 0 . We will take the extremal HS 
held to have a time-independent value on each sublattice, 
h/x and which obey the equations 

hx = ^giXv)^ = 35 ( xa ) ( 6 . 6 ) 


Since the triangles of the kagome can be labelled by 
the sites of a honeycomb lattice on the centers of the 
triangles, we can regard the Hamiltonian of Eq.(6.1) as 


a coupling between the chiralities on a honeycomb lat¬ 
tice. Although Eq. (6.3) has a very complicated form, it 


can be simplihed by using a Hubbard-Stratonovich (HS) 
transformation in terms of a scalar held h{r,t) on the 
sites {r} to the honeycomb sublattice of the triangles of 
the kagome lattice. Upon this transformation, the action 
of the full system, XXZ and chirality couplings, becomes 


S = S 


xxz 


+ J dt-^'^h{r,t)K ^{r,r')h{r',t) 

j, r' 

- J dt'^h{r,t)x{r,t) (6.4) 


where K{r,r') is the coordination (or connectivity) ma¬ 
trix of the honeycomb lattice and K~^{r,r') is its in¬ 
verse. The HS held h{r,t) plays the role of the chirality 
held introduced in Sec |IVB[ except that here it is a func¬ 
tion of time and space. 

We can now apply the hux-attachment transformation 
to a system whose action is given by Eq.(6.4|, and, as 


we did in the preceding sections, map this problem to a 
system of fermions on the kagome lattice coupled to a 
lattice Chern-Simons gauge held. However now they are 
also coupled to the HS helds h{r,t) in the same fashion 
as we coupled the fermions to the chiral operator in Sec. 

WB 


where {x{r,t)) is the expectation value of the chirality 
operator on each sublattice. If we further seek solutions 
that do not break the sublattice symmetry, we obtain the 
simple mean held equation for the chirality 


h = S9 (x) 


(6.7) 


and the critical value of the chirality coupling gc is given 
by the usual mean-held-theory relation 


1 = 35, 


d{x) 


dh 




( 6 . 8 ) 


where 

ah 

model. 


h=0 


is the chirality susceptibility of the XXZ 


For A/J = 0, we hnd that for values of 5 > 13.3J, 
there exist non-vanishing solutions of the chirality pa¬ 
rameter i.e. h > 0 as can be seen in Fig. In these 
cases, we end up with a non-zero chiral term similar to 


that of Eq. (3.1 1 and the resultant phase would again be 


gapped and correspond to the chiral spin liquid discussed 
in the previous section. The critical value of 5 reduces 
as one approaches the isotropic point. For A/J = 1, the 
critical value is much smaller, gc « 3.4 J. Below this crit¬ 
ical value of gc, the value of h that satisfy the mean-field 
consistency equations are h = 0. In this situation, we 
are back to the situation with just the XXZ Heisenberg 
model and the resultant phase at half-filling would be 
gapless. The expectation values displaced in Fig are 
quite small. The main reason for this is that each chiral¬ 
ity operator has a term proportional to S^. When all the 
sites are e xact ly at half-filling this terms is equal to zero 
(See Eq. (2.4)) and the chiral expectation vanishes. In 


We can now integrate out the fermions, we obtain the order to open up a gap, the densities have to be slightly 
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FIG. 9: (Color online) Expectation value of the chirality 
operator (x) plotted as a function of g for A = 0, the 
XY limit (full circles) and A/J = 1, the isotropic 
Heisenberg point (triangles). The mean-field theory 
critical values are gc ~ 13.3 J in the XY limit, A = 0, 
and gc ~ 3.4J at the isotropic Heisenberg point, 
A/J= 1. 


shifted away from zero giving rise to a small non-zero 
chiral expectation value. 

This leaves us with the question of what is the ground 
state of the XXZ Heisenberg antiferromagnet on the 
kagome lattice for small A < 1 and g < gc- Naively, we 
would seem to predict that it is equivalent to a theory 
of two massless Dirac fermions which, on many grounds, 
cannot be the correct answer. In fact, Lopez, Rojo and 
one of ud^ found the same result in the XY regime of the 
quantum Heisenberg antiferromagnet on the square lat¬ 
tice (which is not frustrated). These authors showed that 
the naive expectation is actually wrong and the fermions 
became massive by a process that can be represented 
as the exchange of Chern-Simons gauge bosons. Due 
to the stronger infrared behavior of the Chern-Simons 
gauge fields (compared with , e.g.. Maxwell), this ex¬ 
change term leads to an induced mass term for the Dirac 
fermions which is infrared finite (but linearly divergent 
in the ultraviolet). Most significantly the sign of the in¬ 
duced mass term leads to an extra Chern-Simons term 
which exactly cancelled the term introduced by flux at¬ 
tachment, leaving a parity-invariant Maxwell-type term 
as the leading contribution to the effective action. Fur¬ 
thermore, in 2-1-1 dimensions, a Maxwell term is known 
to be dual to a Goldstone boson. Lopez et al. con¬ 
cluded that the ground state of the antiferromagnet on 
the square lattice in the XY regime has long range order 
and that the Goldstone mode is just the Goldstone mode 
of the broken U(l) symmetry of this anisotropic regime. 


It should be apparent that in our case we can repeat 
the same line of argument almost verbatim which would 
suggest that in the XY regime the ground state of the 
antiferromagnet on the kagome lattice should also have 
long range order with a broken U(l) symmetry. However, 
this conclusion is at variance with the best available nu¬ 
merical evidence which suggests, instead, that the ground 
state is Z 2 spin liquid (of the Toric Code variety). The 
resolution of this issue is an open question. 

In summary, this mean field theory predicts that be¬ 
yond some critical value of the ring-exchange coupling 
constant g, which in this mean-field-theory is typically 
large, the system is in a chiral spin liquid state with 
a spontaneously broken time reversal invariance. How¬ 
ever, below this critical value the mean field theory seem¬ 
ingly predicts that the Heisenberg antiferromagnet on the 
kagome lattice is in a phase with two species of gap¬ 
less Dirac fermions. However this is not (and cannot be) 
the end of the story. Indeed, the fermions are strongly 
coupled to the Chern-Simons gauge field which can (and 
should) change the story. In fact, in Ref. [51] a similar 
result was found even in the case of a square lattice. A 
more careful analysis revealed that, in that case which is 
an unfrustrated system, the fermions acquired a mass in 
such a way that the total effective Chern-Simons gauge 
action vanished, resulting in a more conventional phase 
with a Goldstone mode. At present it is unclear what is 
the fate of the Dirac fermions in the case of the kagome 
lattice. In fact, most numerical data on the kagome anti¬ 
ferromagnet suggests that it is a Z 2 spin liquid. Whether 
a Z 2 spin liquid can be reproduced using our methods is 
an open problem. 


VII. CONCLUSIONS 

In this paper, we investigated the occurrence of chiral 
spin liquid phases in the nearest neighbor XXZ Heisen¬ 
berg Hamiltonian (with and without an external mag¬ 
netic field) on the kagome lattice in the presence of var¬ 
ious perturbations: a) a chirality symmetry breaking 
term, b) Dzyaloshinski-Moriya interaction (only in the 
XY limit), and c) ring-exchange interactions. At the 
mean-field level, we found that in the first two cases 
these interactions open up a gap in the spectrum and 
lead to phases with non-trivial Chern numbers (analo¬ 
gous to an integer quantum Hall state) in the XY limit, 
j < I. When the effects of fluctuations are included, we 
find that these states actually correspond to fractional 
quantum Hall states for bosons with a spin Hall conduc¬ 
tivity of (j'^y = i. This chiral spin liquid state survives 
for larger values of the chirality term but for larger val¬ 
ues of the Dzyaloshinski-Moriya term, the chiral spin liq¬ 
uid state vanishes. Our results qualitatively agree with 
those obtained in a recent numerical study using the same 
mo del .1^ 

We also considered the effects of adding ring-exchange 
term on the bowties of the kagome lattice and found that, 
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provided the coupling constant is larger than a critical 
value (which depends on parameters, e.g. the value of 
the Ising interaction), time-reversal symmetry is spon¬ 
taneously broken and results in a topological state sim¬ 
ilar chiral spin liquid state. However since the critical 
couplings that we find are rather large, ranging from 
^ ~ 13.3 in the XY limit to j ~ 3.4 at the isotropic 
point, we cannot exclude that other phases may also play 
a role. In particular, we have not explored the possible 
existence of topological phases with nematic order.S^ 

In an earlier paper we showed that in the presence of 
a magnetic field, the nearest neighbor Heisenberg Hamil¬ 
tonian gives rise to magnetization plateaus at m = ^, | 
and I in the XY limit.^^ Here, we found that some of 
these plateaus survive with the inclusion of the chiral¬ 
ity and Dzyaloshinski-Moriya terms. In addition we also 
find another plateau at magnetization m = | with a spin 
Hall conductivity cr®^ = |. 

In the absence of an external magnetic field, the flux 
attachment transformation that we use here, at the level 
of mean field theory, naively maps the kagome antiferro- 
magnet onto a system of two species of massless Dirac 
fermions. Since this state is not gapped, the spectrum 
(and even the quantum numbers of the states) is not 
protected by the effects of fluctuations. Of all the fluctu¬ 
ations that are present, only the long range fluctuations 
of the Chern-Simons gauge field are (perturbatively) rel¬ 
evant. Indeed, this problem arises even in the simpler 
problem of the Heisenbergantiferromagnet on the square 
lattice, and Lopez et showed, using a non-trivial 
mapping, that already at the one-loop level the spec¬ 
trum changes from “free” massless Dirac fermions to the 
conventional Neel antiferromagnet with XY anisotropy 
(easy plane). In Sec. |^we derived an effective field the¬ 
ory for the kagome antiferromagnet at zero field and, not 
surprisingly, found a state which naively has two species 
of massless Dirac fermions. A simple minded applica¬ 
tion of the same line of argument would also predict an 
easy-plane antiferromagnet which has a Goldstone mode 
(in the XY regime). This however is not consistent with 
the best numerical data which shows no long range order 
but a topological Z2 state. How to reconcile these two 
scenarios is an open question which we are investigating. 

On the other hand, we should note that, contrary to 
the case of non-relativistic fermions, a theory of massless 


Dirac fermions coupled to a Chern-Simons gauge theory 
is non-trivial. While in a massive phase this coupling 
should also amount to change in statistics, the massless 
case is much less understood. In fact, the only case which 
a related problem is understoocP^"^ is the case in which 
the gauge fields have a gauge group U(N) and the Chern- 
Simons action has level k. In the limit in which N ^ 00 
and fc —>■ 00 (with ^ fixed), this problem maps onto a 
Wilson-Fisher fixed point of a scalar coupled to a Chern- 
Simons gauge theory with gauge group U{k) at level N 
(with the same ratio ^). Away from this regime not 
much is known. In this large N and large k limit, the 
system remains conformally invariant (and hence criti¬ 
cal). Our present understanding of the kagome antifer- 
romagnets suggests that for small enough N the system 
should become gapped and conformal symmetry should 
be spoiled. If the latter scenario is correct, then there 
should be a direct transition from (quite likely) a time- 
reversal invariant Z2 topological phase to a chiral spin 
liquid phase. If this were to hold the quantum phase 
transition would most likely be first order, although an 
exotic Landau-forbidden transition transition is also a 
possibility, perhaps of the deconfined quantum criticality 
type.l^ In the latter case, the above cited recent theories 
of conformal quantum field theories with Chern-Simons 
terms may be natural candidates for the field theory of 
such a quantum critical point 

ACKNOWLEDGMENTS 

We thank Bela Bauer, Andreas Ludwig and Ronny 
Thomale for stimulating discussions. EF thanks the 
KITP (and the Simons Foundation) and its IRONIC 14 
and ENTANCLED15 programs for support and hospital¬ 
ity, and to the Departamento de Ffsica (FCEyN), Univer- 
sidad de Buenos Aires, for its hospitality. This work was 
supported in part by the National Science Foundation by 
the grants PHY-1402971 at the University of Michigan 
(KS), DMR-1408713 at the University of Illinois (EF), 
PHYll-25915 at KITP (EF), and and by the U.S. De¬ 
partment of Energy, Division of Materials Sciences under 
Awards No. DE-FC02-07ER46453 through the Frederick 
Seitz Materials Research Laboratory of the University of 
Illinois at Urbana-Champaign, and Ministry of Science 
and Technology (MINCyT, Argentina). 


^ L. Balents, Nature 464, 199 (2010). 

^ G. Misguich, T. Jolicoeur, and S. M. Girvin, Phys. Rev. 
Lett. 87, 097203 (2001). 

® D. C. Gabra and G. L. Rossini, Phys. Rev. B 69, 184425 
(2004). 

^ S. Nishimoto, N. Shibata, and C. Hotta, Nat. Commun. 
4, 2287 (2013). 

® K. Kumar, K. Sun, and E. Fradkin, Phys. Rev. B 90, 
174409 (2014). 


® S. Yan, D. A. Huse, and S. R. White, Science 332, 1173 

( 2011 ). 

^ H.-G. Jiang, Z. Wang, and L. Balents, Nat Phys 8, 902 

( 2012 ). 

® S. Depenbrock, I. P. McGulloch, and U. Schollwock, Phys. 

Rev. Lett. 109, 067201 (2012) 

“ F. Wang and A. Vishwanath, Phys. Rev. B 74, 174423 
(2006) 

L. Balents, M. P. A. Fisher, and S. M. Girvin, Phys. Rev. 




19 


B 65, 224412 (2002) 

S. V. Isakov, M. B. Hastings, and R. G. Melko, Nature 
Physics 7, 772 (2011). 

G. Misguich, D. Serban, and V. Pasquier, Phys. Rev. Lett. 
89, 137202 (2002). 

G. Evenbly and G. Vidal, Phys. Rev. Lett. 104, 187203 

( 2010 ). 

S.-S. Gong, W. Zhu, L. Balents, and D. N. Sheng, Phys. 
Rev. B 91, 075112 (2015). 

Y.-G. He, D. N. Sheng, and Y. Chen, Phys. Rev. Lett. 

112, 137202 (2014). 

“ Y. Ran, M. Hermele, P. A. Lee, and X.-G. Wen, Phys. 
Rev. Lett. 98, 117205 (2007). 

Y. Iqbal, F. Becca, and D. Poilblanc, Phys. Rev. B 84, 
020407 (2011). 

S. Bieri, L. Messio, B. Bernu, and C. Lhuillier, Phys. Rev. 
B 92, 060407 (2015) 

A. L. Chernyshev and M. E. Zhitomirsky, Phys. Rev. Lett. 

113, 237202 (2014) 

O. Gotze and J. Richter, Phys. Rev. B 91, 104402 (2015) 
M. Freedman, G. Nayak, K. Shtengel, and K. Walker, 
Ann. Phys. 310, 428 (2004). 

A. Y. Kitaev, Annals of Physics 303, 2 (2003). 

E. Fradkin and S. H. Shenker, Phys. Rev. D 19, 3682 
(1979). 

M. P. Zaletel and A. Vishwanath, Phys. Rev. Lett. 114, 
077201 (2015) 

B. Bauer, L. Cincio, B. P. Keller, M. Dolfi, G. Vidal, 
S. Trebst, and A. W. W. Ludwig, Nat. Commun. 5, 5137 
(2014). 

A. Wietek, A. Sterdyniak, and A. M. Lauchli, “Nature of 
chiral spin liquids on the Kagome lattice,” (2015), unpub¬ 
lished, arXiv: 1503.03389 

A. E. B. Nielsen, J. 1. Cirac, and G. Sierra, Phys. Rev. 
Lett. 108, 257206 (2012). 

A. E. B. Nielsen, G. Sierra, and J. I. Cirac, Nat. Commun. 
4, 2864 (2013) 

K. Sun, K. Kumar, and E. Fradkin, “A discretized Chern- 
Simons gauge theory on arbitrary graphs,” (2014), unpub¬ 
lished, arXiv: 1502.00641 
E. Fradkin, Phys. Rev. Lett. 63, 322 (1989). 

D. Eliezer and G. W. Semenoff, Ann. Phys. 217, 66 (1992). 
D. Eliezer and G. W. Semenoff, Phys. Lett. B 286, 118 
(1992). 

V. Kalmeyer and R. B. Laughlin, Phys. Rev. Lett. 59, 2095 
(1987). 


S. C. Zhang, T. H. Hansson, and S. Kivelson, Phys. Rev. 
Lett. 62, 82 (1989). 

J. K. Jain, Phys. Rev. Lett. 63, 199 (1989). 

A. Lopez and E. Fradkin, Phys. Rev. B 44, 5246 (1991). 

K. Yang, L. K. Warman, and S. M. Girvin, Phys. Rev. 
Lett. 70, 2641 (1993). 

X. G. Wen, Adv. Phys. 44, 405 (1995). 

A. Kitaev and J. Preskill, Phys. Rev. Lett. 96, 110404 
(2006). 

M. Levin and X.-G. Wen, Phys. Rev. Lett. 96, 110405 
(2006). 

S. Dong, E. Fradkin, R. G. Leigh, and S. Nowling, JHEP- 
J. High Energy Phys. 05, 016 (2008). 

Y. Zhang, T. Grover, A. Turner, M. Oshikawa, and 

A. Vishwanath, Phys. Rev. B 85, 235151 (2012) 

T. Grover, Y. Zhang, and A. Vishwanath, New Journal of 
Physics 15, 025002 (2013) 

S. Capponi, O. Derzhko, A. Honecker, A. M. Lauchli, and 
J. Richter, Phys. Rev. B 88, 144416 (2013). 

B. K. Clark, J. M. Kinder, E. Neuscamman, G. K.-L. Chan, 
and M. J. Lawler, Phys. Rev. Lett. Ill, 187205 (2013). 

D. J. Thouless, M. Kohmoto, M. P. Nightingale, and 
M. den Nijs, Phys. Rev. Lett. 49, 405 (1982). 

F. D. M. Haldane and D. P. Arovas, Phys. Rev. B 52, 4223 
(1995) 

A. N. Redlich, Phys. Rev. D 29, 2366 (1984). 

E. Fradkin, Field Theories of Condensed Matter Physics, 
Secoond Edition (Cambridge University Press, Cambridge, 
UK, 2013). 

X. G. Wen, F. Wilczek, and A. Zee, Phys. Rev. B 39, 
11413 (1989). 

H. C. Jiang, H. Yao, and L. Balents, Phys. Rev. B 86, 
024424 (2012). 

S.-S. Gong, W. Zhu, and D. N. Sheng, Sci. Rep. 4, 6317 
(2014). 

A. W. Sandvik, Phys. Rev. Lett. 98, 227202 (2007). 

A. Lopez, A. G. Rojo, and E. Fradkin, Phys. Rev. B 49, 
15139 (1994). 

S. Giombi, S. Minwalla, S. Prakash, S. P. Trivedi, and 

S. R. Wadia, Eur. Phys. J. C 72, 2112 (2012). 

O. Aharony, G. Gur-Ari, and R. Yacoby, JHEP-J. High 
Energy Phys. 1212, 028 (2012). 

J. Maldacena and A. Zhiboedov, J. Phys. A: Math. Theor. 
46, 214011 (2013). 

T. Senthil, A. Vishwanath, L. Balents, S. Sachdev, and 
M. P. A. Fisher, Science 303, 1490 (2003). 


34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 



